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Model-independent constraints for the neutron-triton and proton-Helium-3 scattering lengths are 
calculated with a leading-order interaction derived from an effective field theory without explicit 
pions. Using the singlet neutron-proton scattering length, the deuteron, and the triton binding 
energy as input, the predictions a s (t-n) — 9.2 ± 2.6 fm, a t (t-n) — 7.6 ± f.6 fm, a^( 3 He-p) = 
3.6 ± 0.32 fm, and af( 3 He-p) = 3.1 ± 0.23 fm are obtained. 

The calculations employ the resonating group method and include the Coulomb interaction when 
' appropriate. The theoretical uncertainty is assessed via a variation of the regulator parameter of the 

| short-distance interaction from 400 MeV to 1.6 GeV. The phase-shift and scattering-length results 

for the proton-Helium-3 system are consistent with a recent phase shift analysis and with model 
calculations. For neutron-triton, the results for the scattering lengths in both singlet and triplet 
channels are significantly smaller than suggested by R-matrix and partial-wave-analysis extractions 
from data. 

For a better understanding of this discrepancy, the sensitivity of the low-energy four-body scatter- 
ing system to variations in the neutron-neutron and proton-proton two-nucleon scattering lengths 
is calculated. Induced by strong charge-symmetry-breaking contact interactions, this dependence is 
found insignificant. In contrast, a strong correlation between the neutron-triton scattering length 
and the triton binding energy analogous to the Phillips line is found. 
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I. INTRODUCTION 



The refinement of numerical methods (e.g. [Tj-Q) and technological progress made the application of high-precission 
models of the nucleon-nucleon and three-nucleon interaction (e.g. to the four-body system feasible more than 

| a decade ago. By now, the numerical approximations to the solution of the four-body problem are accurate enough 
. to relate the differences in the predicted obscrvablcs solely to the differences in the nuclear force models. Parallel to 
\& ' this development, the effective field theory (EFT) formalism has been successfully applied to an abundance of nuclear 
C***~ . processes, mostly in smaller systems like the deuteron and triton (see e.g. 0, [1[ for reviews). This formalism has 
C*") ■ also been applied in [{J [l(| to understand the Helium-4 nucleus as a universal consequence of peculiar features of the 
I/"") | two- and three-nucleon systems namely unnaturally large scattering lengths and associated shallow bound states. The 
■ results of this work establish for the first time an analogous correlation between the four-nucleon scattering system and 
its subsystems with a nuclear EFT. Furthermore, the absence of bound four-body states in the investigated systems 
' # ■ distinguishes this analysis from previous EFT approaches to the four-body system. 

' Two characteristics of the EFT approach make its application here especially advantageous: first, its ability to 
predict four-nucleon scattering observables with theoretical error bars, and second, the extent to which variations 
in two- and three-body observables flow through into four-body observables can be calculated. Knowledge of both 
the theoretical uncertainty and the dependence on observables in smaller subsystems is required, first, for testing the 
effective range parameters of the four-nucleon system extracted in a phase shift analysis (PSA) on their consistency 
with two- and three-nucleon data. Second, it explains the discrepancy in predictions for those parameters obtained 
from nuclear force models. 

The experimental and theoretical status of low-energy proton-Helium-3 ( 3 Hc-p) scattering is stated briefly to specify 
the focus of this work which is on the singlet/triplet S-wave scattering lengths a^, 4 ( 3 He-p) as a diagnostic for the 

nuclear interaction in this system. The latest experimental extraction of this value via PSA [ll[ yielded: ( 3 Hc-p) = 
11.1 ±0.5 fm and af ( 3 He-p) = 9.05 ±0.15 fm. Both values differ significantly from previous PSAs [HE! (see table H|), 
reflecting the relatively small database especially regarding cross-sectional data below 1 MeV and the associated 
uncertainty in the extrapolation to zero energy. Therefore, naively one would expect that new data could again shift 
the experimental scattering lengths by more than the quoted statistical errors. 

Theoretical predictions are at present inconclusive about the precise values of the four-body scattering lengths, too. 
On the one hand, the model predictions for a^ t ( 3 Hc-p) show a strong sensitivity to the Helium-3 binding energy 

B( 3 He) obtained with the respective model, and on the other hand, a less significant sensitivity to strong isospin- 
symmetry-violating interactions and thereby to the proton- proton scattering length a c (pp). This is manifest in the 
a^ t ( 3 He-p) values of vu and Uig+UIX, where the models differ in _B( 3 He), and in those of V14 and vi$, where the 
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former respects and the latter violates charge symmetry. In table U a selection of recent experimental extractions and 
theoretical predictions is shown. 

To understand this discrepancy amongst theoretical and experimental values the effective field theory approach to 
the nuclear interaction is chosen here. The aforementioned feature of an EFT calculation to provide theoretical error 
bars constrains the observables a^ t ( 3 He-p) to certain intervals. Specifically, the appropriate EFT for very low-energy 
phenomena in nuclear systems, formulated without explicit pions (EFT(jt)), is employed. Its usefulness in describing 
even three- and four-nucleon systems is known by now (see e.g. [8| for a review). With this interaction theory, the 
resonating group method (RGM 3) is used to solve the few-body problem numerically, in which course it is shown 
below that a leading-order (LO) calculation in EFT(^) is consistent with the new PSA extraction for a^ t ( 3 He-p) and 
with high-precision-model predictions for these quantities. 

In addition, the EFT approach is used for a variation of the two-body scattering lengths without losing consistency 
with other low-energy data. A charge-dependent interaction is introduced to shift a c (pp) and the neutron- neutron 
(a s (nn)) scattering length relative to their values obtained in the SU(4) symmetric theory. In that case, a s (nn) is 
equal to the singlet neutron- proton scattering length (a s (np)). It is thus possible to assess the sensitivity of the four- 
nucleon scattering lengths to a s (nri) independently from a change in the triton binding energy B(t) accompanying the 
variation of a s (nn). Regarding experimental uncertainties, the t-n system is more interesting than 3 Hc-p, because the 
experimental knowledge about a s (nn) is comparable to a^ t ( 3 Hc-p), i.e., conflicting extractions stem from insufficient 
data. The proton-proton scattering length, in comparison, is quite accurately determined. It is therefore of interest 
to quantify the effect of a different value for a s (nn) on the prediction for a s / t {t-n) while employing a proper renor- 
malization of the three-nucleon interaction (TNI) so that B(t) stays at the physical value. In this way, the change in 
a s/t {t~ n ) m a ll experimentally admissible scenarios for a s (nn) is explored without losing consistency with other data. 

The results show two similarities to the neutron-deuteron (n-d) system as investigated in fl4| : an underprediction of 
a s / t {t-n) compared to its experimental value, and a correlation of a s (t-n) to the binding energy of the triton (for n-d 
called Phillips line fl5j ) . This work then establishes the four-nucleon scattering system as a universal consequence of 
the two- and three-body system, when so far, only a bound four-nucleon system in the Helium-4 channel was identified 
as an emergent property of the shallow two- and three-body states (see e.g. [IB]). 

This report continues with a summary of the employed nuclear interaction based on EFT(/f) and the specifi- 
cation of the numerical parameters used in the RGM calculation. The following presentation of the results for 
a^, t ( 3 He-p), a s / t (t-n), and the correlation between a s (t-n) the triton binding energy - the analog of the Tjon line fl7j 
for a scattering observable - includes a discussion about the reliability of the error estimates, which is crucial for the 
conclusion. 



II. FORMULATION 



The EFT without explicit pions at LO yields a potential 



l<3 



A\ri-Tj\ 



+ C^f A (f ij ) (1+T il3 ) (l+T 3 - 3 ) 



+C nn (1 - 74,3) (1 - T j>3 ) Mn,)} - (1 - Si ■ Sj) + J2 /a(^) ' httk)Csm n ■ f 3 , 
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cyclic 



(i) 



where a vertex with low-energy constant (LEC) C was regulated via C —> e (p p ) / A C* with p(p') denoting the 
center-of-mass momentum of the incoming (outgoing) nucleon pair ij yielding the radial dependence /a (^3) : = 

( 3,^3/2 ) e~~ Ti i on the relative coordinate fij. The (iso)spin matrices (f)S with indices specifying the nucleon i and 

Cartesian-z-componcnt 3 project onto spin singlet and triplet. The Coulomb interaction with electric charge e repre- 
sents the long-range part of the interaction. For a charge-independent interaction, the three LECs, {Cs, Ct, C3N1}, 
were fitted to a s (np), the deuteron binding energy B(d), and B(t). When charge- independence/symmetry-breaking 
interactions are included, all five LECs {Cg, Ct, C nn , C pp , C3N1} were fitted to {a s (np),a c (pp),a s (nn),B(d),B(t)}. 
The strong isospin-violating terms are suppressed relative to the isoscalar LO interaction and are presumably of the 
same order as the momentum dependent NLO isoscalar vertices [lH ■ For the prediction of a s / t (t-n) and a^ t ( 3 He-p) , 
the inclusion of C nn,pp is another probe of sensitivity to higher-order contributions in addition to the cutoff variation 
and analogous to a replacement of at(np) by B(d) to fit Ct- These two methods are used to map out the LO uncer- 
tainty due to suppressed short-range interactions: a variation of input data and a change of the regulator parameter 
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A. For A, values of 400, 800, and 1600 McV were used. The lower bound was chosen in the vicinity of the natural 
breakdown scale of EFT(/f) while the A dependence up to 1.6 GeV was found sufficient for an extrapolation to A = oo. 

The particular incarnation of the RGM used here was shown to be accurate in its prediction of the four-bod 
scattering system with modern nuclear forces [lj| . For details about the RGM the reader is therefore referred to 
where the method is introduced. Here, only the relevant parameters of the method are listed. The model space used 
to obtain the results presented here includes two two-fragment channels for each four-body calculation, 3 He-p, singlet- 
dcutcron(d )-nn, and t-n, d-pp respectively. In the triplet channels, only the 3 He-p and t-n groupings contribute. For 
LO accuracy in the interaction, it is sufficient to expand the J n = + channel in pure S- waves, i.e., the orbital angular 
momenta on the two fragment-internal Jacobi coordinates of the triton and Helium-3 and on the relative coordinate 
to the neutron or proton are zero. 

The numerical uncertainty due to the finite variational basis was assessed by using two different three-body bases 
(80 and 400 dimensional); by a scaling of the 20 Gaussian widths used to expand the relative wave function by a 
factor of 10; and by including distortion channels for a better approximation of the scattering state in the interaction 
region. The four-body model space was deemed sufficient if the lowest eigenvalue within that space differed by less 
than 60 keV from the respective three-body binding energy. This guarantees that the basis satisfies the minimal 
condition of being large enough to form the bound thrce-nuclcon state with the fourth nucleoli being free. Hence, it is 
expected that the basis is appropriate to expand asymptotic states. For the accurate approximation of intermediate 
states, it was tested if further addition of distortion channels did affect the prediction for the scattering length by 
more than 0.01 fm. Once this level was reached, the basis was assumed to be complete for practical purposes. 

It is worth noting, that the underlying naive expectation that even at LO neither 3 He-p nor t-n sustain a bound 
four-body state is justified by the relatively large energy gap of about 5 MeV from the t-n/ 3 Hc-p threshold to the 
lowest resonant states. In contrast, a similar argument would fail in the Hclium-5 system where it is not obvious 
why the observed shallow P-wave resonances should not become bound at LO in EFT(/f) given the accuracy of a LO 
calculation. 

The variational RGM solution for the S-matrix is obtained for energies below E c . m , = 5 keV and parameterized by a 
phase shift S via S = e 2iS ^ E °- xa ' . To extract a scattering length a s , the standard effective range formula is appropriate 
for t-n, 

(2) 



k cot S(E a . 



while for fragments carrying charges q\.2 its generalization (e.g. pOj ) describes the relative phase shifts in the presense 
of a long-range Coulomb potential 



a 



c 



(Co(r?)fc [cot 6 -i} + 2kr)H(r))) 1 with (3) 



C 2 W = ^ ; ; *= ^755*7; ff (,) = S + J-- Mi „ . (4) 

Formulas (|2l3p were fitted to the respective RGM EFT(^) predictions for the phase shifts where the optimal effective- 
range parameter a^ t ( 3 He-p) / a s / t (t-n) was found stable against inclusions of higher-order terms in the expansions. 

Explicitly, even if — - + §fc 2 was fitted to kcotS the result for a did not shift by more than 0.05 fm. 



III. RESULTS 



First, the effective range parameters calculated with the RGM employing EFT(tt) interactions are shown in tableQ] 
The predictions are compared to previous calculations selected to include models which yield the physical three-body 
binding energy (u 18 +UIX, MT-I/III, and U14+UIX), which only reproduce the two-nucleon sector accurately (t>i4 is), 
and which break (yis) or conserve (um, MT-I/III) isospin symmetry. The given theoretical uncertainties for the 
EFT(/f) predictions of this work are derived from cutoff- and two-body input-data variations. The numerical RGM 
error was assessed independently and contributes a relatively small amount of 0.05 fm to the total uncertainty. For 
the central value, the EFT prediction corresponding to the smallest cutoff was chosen because then the interaction, 
which was renormalized to B(t), yields a Helium-3 binding energy very close to the datum. Increasing A and thereby 
admitting higher-momentum modes in the bound states leads to a stronger contribution from the repulsive Coulomb 
force effectively lowering £?( 3 He). 

The error analysis discriminates between sensitivity of a^, t ( 3 He-p) / a s / t (t-n) to a pp / nn and B( 3 He/t) respectively. 
A shift in the three-body binding energies, induced by a change in the two-nucleon scattering lengths, is absorbed in 
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TABLE I: Theoretical and experimental values for a s / t (t-n) and a^ t ( 3 He-p) scattering lengths. Predictions of charge- 
(as)symmetric interactions with and without three-body models are compared to phase-shift-analysis extractions based on 
different sets of data. 



Theory 



af ( 3 He-p) [fm] 


af ( 3 Hc-p) [fm] 


a s (t-n) [fm] 


a±{t-n) [fm] 


interaction (method) 


Ref. 


11.5 


9.2 


4.10 


3.63 


MT-I/III (FY) 


[25] 


12.7 


- 


4.28 


3.81 


«14 (FY) 


[25] 


11.3 


- 


4.04 


3.60 


i>i 8 +UIX (FY) 


[25] 


11.5 


9.13 


4.06 


3.59 


t)is+UIX (CHH) 


[26] 


12.9 


10.0 


4.28 


3.73 


v ls (CHH) 


[26] 


- 


- 


4.21 


3.54 


t)i4+UIX (CHH) 


[26] 






4.32 


3.79 


vi4 (CHH) 


[26] 


9.2 ± 2.6 


7.6 ± 1.6 


3.6 ± 0.32 


3.1 ± 0.23 


EFT(f) , A — > oo (RGM) 


This work 






Experiment 






of ( 3 Hc-p) [fm] 


af ( 3 Hc-p) [fm] 


a s (t-n) [fm] 


a t (t-n) [fm] 


datasct 


Ref. 


10.8 ± 2.6 


8.1 ± 0.5 






[13] 


[13. 26. 271 


15.1 ± 0.8 


7.9 ± 0.2 






[13,^8,^9], solution 1 


"[12J 


7.2 ± 0.8 


10.4 ± 0.4 






[13. 28. 291 , solution 2 


[12] 


11.1 ± 0.5 


9.05 ± 0.15 






[11. 13. 28-301 


HI] 






4.453 ± 0.1 


3.325 ± 0.016 


[31] 








4.98 ± 0.29 


3.13 ± 0.11 


[32] 





the TNI as the triton is fixed to have its physical binding energy. The effect on a^ t ( 3 He-p) / a s / t (t-n) due to a shift 
in a pp / nn is thereby found insignificant compared to the effect of a change in the binding energy of the three-nucleon 
target. Specifically, the predicted value for a s / t (t-n) changed by less than 1% when C nn is fitted to either of the two 
proposed experimental values for the neutron-neutron scattering length, i.e., a s (nn) = —18.7 ±0.7 fm [2ll |22| and 
o s (jm) = -16.1 ± 0.4 fm [H as long as the TNI is used to fix B{t) = 8.48 MeV. The model predictions are in 
agreement with this finding of a stronger dependence of four-nucleon scattering parameters on the target's binding 
energy than on two-body scattering lengths. The EFT values for a^( 3 He-p) and ap( 3 He-p) are consistent within 
error bars with model calculations which yield the physical three-nucleon binding energies. While the ~EFT(jt) results 
for 3 He-p agree also with the latest experimental numbers, in the t-n system, the EFT(jt) predictions are significantly 
smaller than suggested by models and R-matrix data. The discrepancy to the models is explained below with the 
correlation between B(t) and a s / t (t-n) and thereby with significantly different three-body binding energies. It is also 
interesting to note that for the difference a s (t-n) — a t (t-n), all theoretical calculations, including this work, predict 
approximately the same number and hence collectively differ from the experimental extractions. 

In fig. [TJ the above mentioned correlation between B(t) and a s (t-n) is shown to resolve the discrepancies in the 
theoretical predictions. For each cutoff, the TNI parameter is varied generating the correlation lines (solid black for 
A = 400 MeV). The lines map out a band (gray) whose width at fixed B(t) represents the LO uncertainty of the 
calculation. Again, A = 400 MeV is chosen to generate the central line. The model predictions (empty symbols) 
are consistent with this analysis as they arc included in the band. The correlation explains the difference in model 
predictions for the four-body scattering length as being largely due to their difference in B(t). In the sense that a 
four-body observable is correlated to a three-body observable, the finding is analogous to the Tjon correlation [ljj 
between the two binding energies B(t) and £?( 4 Hc). By correlating a scattering to a bound-state observable, the graph 
shown in fig. Q] is the analog of the Phillips line (empirical: [l5j, EFT explanation: [33l l34l]). 

It is worth mentioning as an outlook that, like the Tjon line, the above correlation is support for a conjecture 
formulated in [3oT ] that the parameters of quantum chromo dynamics place it close to a critical renormalization-group 
trajectory. A small variation in those parameters, namely the up and down-quark masses, would then yield divergent 
two-nucleon scattering lengths and a universal Efimov spectrum in the three-nucleon system. In this socalled resonant 
limit, the analogous bosonic system was shown both theoretically [HI, [37j and experimentally [III H^] to have two 
bound four-body states associated with each Efimov trimer. The bound state and a shallow resonance in the Helium- 
4 + channel hint towards an analogous universal four-fermion system emerging from three resonantly interacting 
nucleons. Whether or not a bound Hydrogen-4 and/or Lithium-4 nucleus will emerge in the limit of diverging two- 
nucleon scattering lengths would therefore further the understanding of atomic and nuclear systems within the same 
theoretical framework. 

Finally, the predictions of EFT(^) for the singlet and triplet S-wave phase shifts for the 3 He-p system are compared 
to the PSA of [ll| in fig. [5] The LO results of this work constrain the phases to the plotted band which was mapped 
out by the cutoff- and input-data variation. The upper edge of the band corresponds to A = 1.6 GcV while the band 
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B(t) [MeV] 



FIG. 1: The singlet S-wave neutron triton scattering length as a function of the triton binding energy. The correlations for 
EFT( f) with regulators from 400 MeV to 1.6 GeV are compared to R-matrix data [U (red solid lines), and to the MT-I/III [U, 
V14 [S]], and the vis+UIX [f| force models with a 3 (t-n) from [25l |. 




012345 012345 
E c . m . [MeV] B c .». [MeV] 

FIG. 2: Singlet (left) and triplet (right) S-wave phase shifts for 3 He-p scattering. The LO EFT(tQ error band resembles a 
cutoff variation between 400 MeV (solid line, central phase) and 1.6 GeV (upper edge of shaded band). Included are the most 
recent PSA values [ll[ (points). 

center is defined with A ~ 400 MeV. The width of the band increases with energy, consistent with the expectation of 
a larger uncertainty due to a slower convergence of the EFT at higher energies. Note that if the lower edge is used 
as the central phase shift, the PSA values lie within the uncertainty range. No significant shift in the phase shifts 
over the considered energy range was observed when the proton-proton scattering length was fitted by considering 
charge-symmetry-breaking interactions. The spread is induced predominantly by the varying _B( 3 He) as only B(t) is 
held fixed. 



IV. SUMMARY 

Model-independent constraints for the singlet and triplet S-wave proton-Hclium-3 and ncutron-triton scattering 
lengths were calculated. The strong dependence of the predictions on the binding energy of the respective three-nucleon 
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target was quantified by a calculation of a leading-order correlation between the triton binding energy and a s / t (t-n). 
On the other hand, the dependence of the four-nucleon scattering lengths on their two-nucleon siblings a s (nn) and 
a c (pp) was found relatively weak. The effcctive-field-theory approach ensured consistency of the predictions, first, 
with relevant two- and three-body data, and second, with the most recent phase-shift analysis in the 3 He-p system. 
A puzzling underprediction of parameters remains in the t-n case, which is also found when sophisticated models are 
used for the prediction of t-n scattering lengths. 
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